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Abstract 

In this work we calculate the four-point correlation function of vector quark currents 
of QCD via holographic QCD model. Computing the correlator we take into account 
the exchange of vector and axial vector bosons and dilaton in the bulk. The result is 
used for calculation of the two-point correlator of electromagnetic currents in external 
magnetic field at zero momentum, related to RHIC experiments, chiral magnetic effect 
and lattice study. At zero temperature we find this quantity to be loosely connected 
with chiral symmetry breaking and strongly dependent on the confinement properties. 
Some features of the AdS/QCD models are also discussed. 

1 Introduction 

Nowadays the holographic approach to QCD is a common tool for exploring the strong 
interaction properties at nonperturbative region. There is a number of models [TJ 121 El IH , 
which can reproduce meson decay constants, Regge trajectories of meson masses, phenomena 
of spontaneous symmetry breaking and chiral anomaly, confinement of quarks etc. At the 
same time the model which would incorporate all these features simultaneously, does not 
exist. Although holographic models can not yet provide a precise calculation tool, they 
can give some new hints about the nature of the considered phenomena and reveal the 
connection between different features of the theory. Moreover, the features of the model 
should be studied to clarify the further modifications. 

The paper concerns the calculation of the 4-point correlator of vector currents of QCD 
in the framework of the holographic model containing the vector, axial vector, scalar and 
dilaton dynamical fields. Similar strong coupling calculations were made in different models, 
where the photon correlators in supersymmetric Yang-Mills or conformal QCD [HI UJ |8] and 
various features of low energy QCD j9l UHl [TT] were studied. We assume that scalar and 
dilaton field have some arbitrary background profiles, moreover we incorporate the arbitrary 
deformation function of the metric in the infrared region. Thus we get the result, applicable 
to a large class of particular holographic models which appear in the literature. For the 
perturbative calculation we adopt a convenient Witten diagram technique to represent the 
action on the classical solutions as the sum of tree diagrams. 
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The second goal of this work is to study the electric current magnetic susceptibility, 
namely the electric current two-point function in the external magnetic field at vanishing 
momentum. This entity may be related to the measured charge asymmetry in heavy ion 
collisions at RHIC experiment [12], and is studied extensively According to [12], the suscep- 
tibility is generated in the moment of collision by the chiral magnetic effect and the chiral 
chemical potential. The other input is produced by the nonzero temperature of the fireball. 
We study this object at zero temperature and zero chiral chemical potential. Though this 
setting is not sufficient to study the phenomena in the ion collision, we try to find out if any 
other origins of current susceptibility are present in nonperturbative QCD. We see that the 
susceptibility in the holographic model is nonzero at quadratic order in the field, and it is 
directly related to the 4-point correlator computed in the first part. 

As the main result we make some conclusions about the behaviour of the quantity of 
interest. It turns out that the susceptibility at zero temperature behaves as Q A at small 
Q. It is not very sensitive to the pattern of chiral symmetry breaking present in the model, 
namely it is nonzero even when the quark mass and condensate are taken as zero, but depends 
strongly on the way the conformal symmetry is violated, or of the confinement properties. 

As to the perturbation theory in the bulk space of AdS/QCD, we find that usually the 
exchange of the dilaton makes no difference being suppressed by the factor of N c . Also we 
find that the exponential growth of the bulk-to-bulk propagators in the soft wall model with 
positive dilaton makes the perturbation theory there inconsistent. 

In Section 2 we present the setting we are working in, describe the principles of the 
diagram technique and calculate the four-point function of the electromagnetic currents. 
In Section 3 we use this result to compute the electric current susceptibility in various 
models and investigate its dependence on the parameters. The last section is devoted to the 
discussion of the results obtained. 

2 The vector current four-point function 

2.1 The setting 

The AdS/QCD model used in this work contains vector (V®) and axial vector (A*) fields 
in the adjoint representation of the gauge group, which is dual to vector currents of QCD 
((qt a j^q) and (q^^j^q) , respectively), the bifundamental scalar field (X alS ) representing 
the scalar quark current {q a q^), and a dilaton, which corresponds to the trace of the gluon 
field strength tensor squared ( Tr(G 2 ) ). The gauge group is the flavor group U(Nf) of QCD 
and the model is non-Abelian. Generators of U(Nf) are normalized the standard way: 

[t a ,t b ] = if abc t c , (t a t b ) = ^5 ab . 

The model is built in 5-dimensional space, which has an AdS metric with equal to unity 
curvature radius and some arbitrary warp factor W(z). It can be cut at some infrared 
brane, located at z — z m . W(z) tends to zero at the boundary, in order the space to be the 
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asymptotically AdS. 



e 



-2W(z) 



d s = — (— dz + dx^dXfj,), z < z Tl 



The action is a sum of Yang-Mills and Chern-Simons parts. We study the dynamics of 
the fluctuations of scalar field around its vacuum profile X (z) = \x( z ), represented by 
the pseudoscalar field n, X = X^e™ . The fact that only axial vector field interacts with 
this background is a sign of the chiral symmetry breaking in the model. The dilaton field 
can have a vacuum profile as well. We will denote it by &{z) and consider the fluctuations 
4>{Qi z). In the following it is convenient to introduce the sum of the warp factor and dilaton 
vacuum profile: &(z) = W(z) + &{z). The pure Yang-Mills part is symmetrical with respect 
to the exchange of left and right currents, the CS part is anti symmetrical. 

S = J fxdzy/ge-™ (-^) <i$ + F\ + 6 X (z) 2 (A - Ott) 2 + A 2 ^) 2 + M(0)> (1) 

+ (Jt^) (VAF V AF A + VAF A AF V + AAF V AF V ) 

Here the trace is taken over the gauge group U {Nf), A is a normalization constant of dilaton 
field, which should be fixed. M(<f) is a possible potential for dilaton, which is not essential in 
our calculation, as long as we are not interested in the gluon condensate corrections. All the 
other constants in the action are fixed by matching with QCD sum rules (see [U EJ [131 ED- 
In the vicinity of the boundary of the AdS the asymptotics of the scalar vacuum profile is 
x{z) = mz + az 3 . The normalization is chosen so that m is the quark mass in QCD. In the 
action we have omitted all terms with interactions unnecessary for our calculation. These 
are interactions of the second order in axial field in the Yang-Mills part and terms of fourth 
and higher order in fields in the Chern-Simons part. 

Before we proceed, we have to fix the constant A in the action. This can be done by 
matching with QCD sum rule result for the gluon correlator 



j d*x(G 2 (0)G 2 (x))e^ = -^LAqHu (Opj . (2) 

In order to compute this correlator in the holographic model it is enough to solve the lin- 
earized equation of motion for the dilaton near the boundary of AdS 



d 2 <p - - z d z 4) + Q 2 4> = o. 



The solution (bulk-to-boundary propagator) is 

^\z) = ^z 2 K 2 {Qz) + C 2 z 2 I 2 (Qz)}, 

where / and K are the modified Bessel functions, and is a boundary value of dilaton field, 
which should be associated with the source of corresponding operator in the generating 
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functional of 4D theory, namely Tr(G 2 ). The final step is to insert this solution into the 
action and take the second variation with respect to the source 4>{Q). We get the two-point 
function at large momenta Q: 



d4>(Q)d^(-Q) V24vrV 8 

Comparing this with sum rule (T5]) we see: 



,2 48(JV C 2 - 1) 



* = (3) 



' C 



2.2 The variation of classical action 

We proceed to our main task, the calculation of four-point function of vector currents. 
According with the general AdS/CFT recipe [16] we take the variation of classical action 
with respect to four boundary values of vector fields 



(JZ(h)J b a(k2)J'(h)Jt(h)) = —~ ; » * S cl (V 4 



So we need to calculate the classical action up to the fourth order in boundary val- 
ues, which involves the calculation of classical solutions up to the third order in sources. 
The calculation can be done perturbatively. Let's make a Fourier transformation over four 
coordinates and rewrite the action ([1]) in a more explicit form: 

s = s v + S* + S A 



'24V 

~n7 



S v = I d%d% S\q 1+ q 2 ) —(-)V^ qi ,z)d z V^q 2 ,Z 

i 



1 .2=0 



+ / d%d%dz 5\ qi +q 2 ) V2(qx,z) 



z 



z „ e 



^d z ^d z g atl - q 2 g a/1 + q 2a q 2fl 



+ / d 4 qi d%d%dz (2i)T a XV-( qi ,z)V^q 2 ,z)V 7 c (q 3 ,z) 



+ / d%d%d%d%dz — (-)^Q^ V:( qi ,z)V^q 2 ,z)V^q 3 ,z)V s d ( qi ,z) 



(4) 

(5) 
(6) 

(7) 



In this term the partial integration has been used in the kinetic term, and we use the notation 
for the triple and quad vertex functionals: 



T $7 = ^(qi+q 2 +q 3 ) / a6c (#a/39l 7 + 9fh(ha + 9^Qw)i 



vribcd ;4 



S\qi+q2+q3+q4) ^^9^905 



(8) 
(9) 
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Similarly, the partial integration is made in the dilaton part of the action, but due to the 
absence of the dilaton source in the problem, the boundary term vanishes, and only kinetic 
and interaction terms survive 



24vr 5 



,4 ,4 , 48(iV c 2 -l) e" 
a qid q 2 dz — (p{qi,z) — 



-<S>-2W 



-&-2W 



-&-2W 



-d z - q\ 



+ / d A q\d A q 2 d A q 3 dz 



l b ^(q 3 ,z)V:( qi ,z)V»(q 2 ,z). 



(10) 

(11) 



The notation for dilaton vertex functional is 



D 1 



ah 



S\qi+q2+q 3 ) { - S ab dld 2 z g a ^ + S ab (-( qi q 2 )g^ + q w q 2a )}, (12) 



where by definition momenta (91,52,93,94) and differentials dl,d 2 , . . . correspond to fields 
with flavor indices (a, b, c, d) respectively. 

The CS term warrants more careful study. Its explicit form in the coordinate space is 



/N 
d " XdZ Y^ Mt a tH c )e^(d z A^[d a V u a V p b - V u a d a V b ] 

+ Al[d z V v a d a V b - d a V y a d z V b ] 
+ 9^[W;-^]). 

We perform the integration by parts in z coordinate and we omit the boundary value, because 
the axial source is not present. So in the momentum space the part S A looks like 



24?T' 



d 4 q 1 d A q 2 dz 5 A (qi+q 2 )x 



x \ K{<h,*) 

-&-2W 



- d z d z g afl (q 2 g af + - 92a92 M ) 



A a Jq2,z) 



(13) 



+ 3- 



-X(z) 2 z) - q lfl n( qi , z)) (A«(q 2 , z) - q 2f ,7v(q 2 , z)) + (d z n(q, z)f 



+ J d A qi d A q 2 d%dz (24*) A^ y V:( Zl , qi )V b (z 2 ,q 2 )A°(z 3 ,q 3 ) 
where we've introduced the notation for the axial vertex functional 



* = 5 4 (9i+92+93) (t a t b f)e a ^(d l z q 2 a - d 2 z q l a ). 



(14) 



(15) 



Now we should insert classical solutions calculated order-by-order in vector field sources: 
V = + V® + V^; A = A®] <f> = (f) {2) . One can check, that the term with V^V® in © 
cancels the boundary term (j4"l). So we are left with the four different diagrams of Fig.l. (la) 
comes from the connection of two via kinetic term in (jSJ) and connection of two 
and one by the triple vertex in ([6]). The connection of two A^> in the quadratic term of 
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Figure 1: Tree diagrams, corresponding to the classical action of the fourth order in vector 
sources: a) vector intermediate boson, b) axial intermediate boson, c) exchange of dilaton, 
d) quad vector vertex 



(fT5D and the vertex of jUD with two V w and one contribute to (lb). Similarly (lc) is 
made of two ^ in ([IU]) and two V {1) and (2) in (HI]). The dia gram (Id) is the connection 
of four via the quad vertex in (J7J). 

We proceed to the calculation of classical solutions. Let's start with the vector field 
V^(Q, z). Its equation of motion is 



d z d z g afl H (<?i5V - giaft/J 

z z 



V a A-qi,z) 



d%d%iT a a %V(q 2 ,zW c (q 3 ,z) 



As we are interested in and only, we've dropped unnecessary interaction terms. 
The solution, which is linear in source, is obtained by means of bulk-to-boundary propaga- 
tor v(Q,z), which is the solution to the linearized equation of motion with unit boundary 
condition at z — > 0, v(Q,0) = 1. Using this, we write out the expression for the classical 
solution: 



V^ a (q,z) = V:(q)v(q,z) 



(16) 



The second order solution can be found perturbatively with the use of the Green function 
of the equation. Due to the nontrivial vector structure of the equation, we introduce two 
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Green functions, the transverse G ± (Q, z) and the longitudinal one z) with respect to 

the momentum Q. They are solutions to the following equations: 



d z — d z + — q\ 

z z 



G ± (q u z,z') = S(z,z'), 



(17) 



d z —d z G^ qi ,z,z')=S(z,z'), 



and govern the iteration procedure. For short we will denote their sum by G\ x , = [gw — 
® X q 2 x ' ] G ± + ® X ® 2 X ' GK Now we can straightforwardly find the second order solution V^: 

V<® a a (-<h,z) = [ d%d% iTf^V^V^qs) x / dz'^Gl al ( qi ,z,z')v( qi ,z')v(q 2 ,z'). 



The equation of motion of dilaton field is 



d z 



-d z + 



e 2 



<p(-qa,z) 



1 e 



^V:( qi ,z)V^ Q2 ,z). 



2A 2 



As only sources for vector fields are present, we calculate the second order solution immedi- 
ately. We use the Green function defined by the equation 



~<S>~2W 



o. 



-3*4- 



$_2vy 

e 2 



G*( q3 ,z,z>) = 5(z,z>) 



and get the solution: 



d qi d q 2 



2A 2 



dz'—G*(q 3 , z, z')B%v( qi , z')v(q 2 , z')V:( qi )V>( q2 ). (19) 



The last part is to compute the solution for axial vector field from the equation 



_cj> _<T> 

e „ e 



d z d z g lfl + (g 3 ^ 7M - g 37 93/J 



-3>-2W 



+ 6- 



X {z) 2 (A*{- qz ,z) + q^i-qs, zj) = 12z j d%d% K^V^ qi )V b p {z,q 2 ). 



(20) 



This equation is more complicated than in the vector case, due to the interaction with a 
scalar profile. We will compute the explicit form of the Green function later. Nevertheless 
we can use the same notation for its longitudinal and transverse parts as for the vector field 
and denote it Gw. The solution is 



77 



A%(-q 3 ,z) = / dz'd%d% (12i) A^ 7 ,V; a ( gi )Vj(g 2 ) G^,(z, z' , q 3 )v(z' , qi )v{J, q 2 ) (21) 
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Having the solutions (fT8"l) . (ll9p and (I2T1) at hand, we are able to write out the classical 
action of the fourth order in sources. It is 

= j d A qid%d%d%d% 8\ qi +q 2 +q 5 ) 5\q 5 -q 3 -q,) V^V^V^Vg^X 



(-1)-^ T%T^ 5 x / dzdz' r GUq 5 ,z,z')v( qi ,z)v(q 2 ,z)v( q3 ,z')v(q,,z') 



Nc f , , .. e-*e"< 

4 2 



(22) 

4 24vr 2 48(iV 2 -T) ) D S D 5 x V V^^ 5 ' *' *>(^M^ ^Mfe.^Cfli.^) 

(23) 



_ 6 _| y ^ A^A^ x Gi xl (z,z',q 3 )v(z', qi )v(z',q 2 )v(z,q 3 )v(z,q 4 ) 



(24) 



>W x J dz—v(qi,z)v(q 2 ,z)v(q 3 ,z)v(q 4 ,z)j, (25) 

where the following structures appear (we omit here delta-functions, moved to the action): 

K%^s = f abe f ecd (g a pqu + gpeQsc + 9^p) (g^Qss + g-ysQse' + gesq^); (26) 

DjD$ = 5 a6 5 cd {9 2i a Z2 ^ + (gifc)^ - gi^g 2a }{^ 3 ^ 4 ^ + (<m 4 )<? 7 5 - qm^}\ (27) 

A? A A^ A , = (t a tY)(t c t d t e )e^e^(dlql - - (28) 

QX = /^.TW/m- (29) 

Now taking the variation with respect to four sources resolves into the summation of all 24 
permutations of indexes and corresponding momenta in the above structures. Taking the 
variation with respect to the transverse source will result in the transverse structure in the 
correlator. We will not write these projectors expicitly and will imply, that the result is 
transverse with respect to momenta of the corresponding external sources. There is no point 
if doing the variation in general case, so we will make it for the particular application. 



3 Electric current susceptibility 

3.1 Kinematics 

One of the interesting problems is the calculation of the two-point function of vector currents 
in the external magnetic field at vanishing momentum Q. In the first order of magnetic field 
it is related to the correlator of three vector currents, which is zero in our model. To the 
second order in the magnetic field it can be expressed as 

(MQ)M-Q))b = (MQ)Jf}{-Q)Mki)Js(fo))effl klM -M, (30) 

where el, are polarizations of the real photons. We encounter the 4-point correlator, which 
we can now calculate. As it involves only electromagnetic currents ( J em = | J° + | J 1 ) two 
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terms in the variation vanish immediately. These are (|22|) and ( )25|) . because they contain 
structure constants f abc , which are zero if one of the indices is zero or if any two are identical. 
Moreover, we see, that the term, that contains intermediate dilaton (123"j) is suppressed by 
the factor of N c , comparing to the others. So it is also negligible. 

The main input to the correlator is from the axial boson interchange fl24|) . Let's take 
the sum over all permutations of sources V£ (Q), Vf(— Q), V^(ki), V d {k 2 ) in the structure 
(1281) . Due to the symmetry properties of this structure all 24 terms form 3 groups of 8, 
corresponding to S-,T- and U- channels: 

S-channel, 2({t a t b }t e ) ({t c t d }t e )e a ^e^ x '^(d^Q 2a - d^Q la )(d k Jk 2n - d k Jk l7T ) 
T-channel, 2({t a t c }t e )({t b t d }t e )e a ^e^ x '^(df 1 k 1(7 - d^Q la ){dfk^ - d k JQ 2w ) 
U-channel, 2({t a t d }t e )({t b t c }t e )e a5x °e^ x '*(d^k 2a - d k /Q la )(d k JQ 2 „ - dfk ln ). 

Here the superscript of the differential df denotes the particular bulk-to-boundary propa- 
gator v(Q,z) that should be differentiated. In these expressions we keep only terms linear 
in k\, k 2 to study the dependence on the magnetic field. Consequently, the S-channel contri- 
bution vanishes. We multiply the result by the polarizations of real photons and figure out 
the dual strength tensor of external field e^ upT k p e T = ^F^ u . Finally, we get 

T-channel, 2({t a t c }t e )({t b t d }t e ) Qf *^) F^ x '(k 2 )j d^d^\ 
U-channel, 2({t a t d }t e }({t b t c }t e ) (-F aX (k 2 ) F? x ' (h) ) 3^d Q / . 



Our main task is to compute the correlator of electromagnetic currents ( J em = | J° + | J 1 ] 
We can express it via our obtained four-point function: 



(J em {Q)J em {Q)J em {k)J em {k)) =—(J J°J°J ) 

1 1 



+ 2~[(J 1 J°J J ) + (J J°J 1 J )} 

A 6 

+ 2~[(J°J 1 J 1 J 1 ) + (J\J\J J 1 )] 



6 4 



[J 1 J 1 J 1 J 1 ) 



Having this structure, we can sum up U- and T- channels and get the expression for two-point 
function of electromagnetic current in the magnetic field. 



(r m {Q)r m {Q)) B 



1 111 

h 6 1 

2 4 T 22 g2 T g 4 



l -F aX {k 2 ) F^'ih) 



(31) 



7T- 



dzdz' Gi x ,(z, z', Q)d z ,v(z, Q)d z v{z, Q).v(z, 0)v(z', 0). 
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The number in square brackets comes from the charges calculation and hereafter will be 
denoted by C. 



3.2 The calculation of the coefficient 

To proceed, we have to specify the model we are working in. The simplest choice is a "hard- 
wall" model proposed in [TJ Other possible choices will be discussed further. In the 
"hard-wall" model there isn't any background dilaton profile or warp factor in the metric, 
so we put $ as zero. The chiral symmetry breaking is described by the scalar profile 

x(z) = mz + az 3 

in the whole bulk space up to the IR boundary z m . In the limit of Q — > the equations of 
motion for transverse and longitudinal parts of axial vector field coincide, so the bulk-to-bulk 
propagator can be expressed as 

Gf x ,(Q,z)=g xx ,G{z) + 0{Q), 

where G(z) is defined by the equation 

1 „ 1 



d z -d z - Q^x(z) 2 
z z 6 



G(z,z') = 5(z,z'), 



with usual hard-wall boundary conditions: 

G(0,z) = 0, d' z G(z,z%'= Zm = 0. 
The vector bulk-to-boundary propagator can be found explicitly, it's revealed as: 

v(Q, z) = ^ Zm ] Qz Xi(Qz) + QzJd(Qz), 

where X and /C are modified Bessel functions. One can check, that at Q — > limit v(Q, z) = 1 
and d z v(Q, z) ~ Q 2 , so in the integral fl3Tj) only differentials of propagators remain, which 
result in the overall factor of Q A . Now we can calculate the integral numerically. For the 
parameters z m = (323Mew) _1 , m = lOMev, o = (460Mev) 3 (see appendix of [TJj) we 
obtain: 



[J' m (Q)J' m (Q)) B =C 7 F" A (fe) F*(h) x (-6)^Q 4 (-1.22Ge !) - 



w ' V^° i(l ' 2) P ' ih) ) X ( " 6> 5 

The value of the integral here is computed numerically, it is dimensionful and warrants some 
more study. In the hard wall model we have 3 dimensionful parameters: m,a and z m . One 
can check, that the quark mass has very little influence on the result, so we can freely take 
the chiral limit (m = 0) with the change of about 1%. The a plays a bigger part, but it is 
also not crucial for the result. Taking it equal to zero gives the value 2.93Gev~ 6 and larger 
a decreases this number until at very large a ~ lOGeV 3 the integral practically vanishes. 
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Important for us is that the current correlator in the external field is nonzero even in the 
chiral limit and at zero chiral condensate. This indicates, that its origin is not the chiral 
symmetry breaking. In the case of m = a = the result obviously behaves as z^, or Aq CD , 
so taking the scale parameter to infinity kills the result as well. But concerning the limit 
Aqcd — > we do not find such a smooth behaviour as in the case of condensate. In this 
limit the integral diverges strongly, which is a sign of qualitative dependence on the Aqcd- 
This lets us conclude, that the most important for the considered value feature of the model 
is the IR cut-off the AdS space, related to the confinement or conformal symmetry breaking. 

According to this issue it is interesting to study the "soft-wall" model, because it has an 
improved behaviour in the IR which includes the proper description of the meson masses 
[31 15]. Moreover, we are not very interested in chiral symmetry braking, which is problematic 
in considered framework (see [1]). In this model one has an exponential factor e cz in the 
action, due to the nontrivial dilaton profile $ [5] or to the deformation of the metric W [5]. 
In our case this makes no difference, because it is their sum $ which appears in equations 
of motion. The most popular version of the "soft-wall" model has negative c OH], so that 
the action has a rapidly falling prefactor e -7 z . This gives rise to the strong problem in 
our calculation, because the bulk-to-bulk propagator, being the inverse of the kinetic term 

2 2 

in the action, acquires a factor e 7 * and diverges extremely fast at large z. This makes the 
diagram under consideration 021]) diverge too, because the vertices coming from the Chern- 
Simons term do not have any exponential suppression, due to the independence of metric 
and dilaton. As a consequence, one can not use the perturbation theory here. This problem 
is not surprising, since the model with negative c does not express the confinement (namely 
the proper heavy quark potential) crucial for our computation [T7] ITS] . 

In order to incorporate the confinement, one has to use the model with positive c [T7] [TS] 
[19] . Then all the solutions to equations of motion decay in IR as well as the bulk-to-bulk 
propagator. So the perturbation calculation converges and we can get the striking result: 



It is nonzero at Q — 0, what seems to be quite unphysical in the framework of zero temper- 
ature and absence of any chemical potentials. The reason is the presence of the unphysical 
massless mode in the vector sector in the model with positive c, which has been observed in 
[5] and discussed in [Tj5]. This feature prevents us from trusting the answer above, until the 
nature of this mode is clarified. 

4 Conclusion 

In this work we've performed a perturbative 5D calculation of the vector current four-point 
function of QCD via the holographic model. The generality of the used setting allows us to 
study some common features of such calculations in various AdS/QCD models. Namely we 
find, that due to the normalisation of dilaton field, the exchange of dilaton in the bulk is 
suppressed by the factor of N c . Even though the form of the dilaton field and its features 



(J em (Q)J em (Q))B 



sw 




\F a \k 2 ) Ft 



(h) x (-6)^(0.7G , e^)- 2 
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in holography are still subjects of investigation, we claim, that this behaviour has to remain 
in all the models. The obtained result is applicable to a large class of commonly used 
holographic models and can be used in various problems, as it includes the chiral symmetry 
breaking term, the anomaly term, the arbitrary number of quark flavors, the arbitrary dilaton 
and scalar potentials and the arbitrary warp factor in the metric. 

We use the vector current four-point function to calculate the two-point correlator of 
electromagnetic currents in the external magnetic field. This quantity turns out to be weakly 
dependent on the chiral symmetry breaking in the model and is sensitive to its infrared 
behaviour, namely to the way the conformal symmetry is broken. Studying this calculation 
in the hard wall AdS/QCD model and soft wall models with different dilaton behaviour, we 
find that the perturbation theory is divergent in the soft wall model with falling exponent in 
the action. This problem is not visible until every term of the action, including interactions, 
is multiplied by the same factor. But considering the Chern-Simons action, which is free of 
metric as well as dilaton input, we get the exponentially growing bulk-to-bulk propagator, 
not suppressed by the interaction vertices and leading to the strong divergence of the result. 
This turns to a handicap for calculating diagrams with bulk-to-bulk propagators and Chern- 
Simons vertices in such a framework. In the model with growing exponent in the action 
we obtain the value, finite in the limit of small Q, which is unphysical at zero temperature. 
This results from the presence of the massless pole in the vector two-point function, which 
leads to the violation of charge conservation and is forbidden in QCD. We hope, that the 
nature of this pole and the solution of this defect will allow more accurate calculations in 
this model. 

In the end of the day we obtain the best result in the framework of the "hard-wall" 
model. The value of the correlator in the magnetic field is proportional to B 2 at small 
B, where B is the field strength. At small Q it goes to zero as Q 4 , what is consistent 
with general expectations at zero temperature and chemical potentials. This confirms that 
the holographic model, featuring chiral symmetry breaking and chiral anomaly, does not 
autonomously develop the chiral chemical potential, needed for the nonzero result. Some 
other components are obviously necessary to describe the phenomena of heavy ion collisions. 
Constructing a holographic model which would feature nonzero temperature and topological 
charge of QCD could shed light on this problem. 
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